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Abstract 

We discuss the reported detection by the MILAGRO experiment of localised hot 
spots in the cosmic ray arrival distribution and the difficulty of interpreting these 
observations. A model based on secondary neutron production in the heliotail is 
shown to fail. An alternative model based on loss-cone leakage through a magnetic 
trap from a local source region is proposed. 
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1 Introduction 

The MILAGRO experiment pQ has reported the surprising detection of two 
localised "hot spots" in the cosmic ray arrival distribution the stronger and 
more localised of which (region A) coincides with the local heliotail direc- 
tion. While the intensity is quite low (fractional excess of order 6 x 1CT 4 in 
region A) the observation of what appear to be essentially well-collimated 
beams of hadrons (with energy about 10 TeV or more) coming from very 
specific directions is hard to accommodate within standard models of cosmic 
ray propagation. Although the observed tiny small-scale anisotropy could be 
most naturally explained as being due to gamma-rays, the MILAGRO group 
claims that an electromagnetic origin of the signal is excluded with high con- 
fidence. This leaves protons and nuclei, and, to a lesser extent, neutrons as 
particles which could cause the observed anisotropy. However as noted in the 
MILAGRO discovery paper protons at this energy only have gyroradii of order 
0.005 pc m 10 3 AU in the estimated local interstellar magnetic field of 2 fiG 
and the decay length of neutrons is only 0.1 pc. These simple, but robust, esti- 
mates indicate that any interpretation of this result will require nonstandard 
ideas and approaches. 

That there should be some structure in the arrival direction distribution is very 
natural. Apart from large scale gradients related to the bulk Galactic prop- 
agation, it would be very strange if there were not some signature of nearby 
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recent cosmic ray source events (such as local supernovae) as discussed by, 
for example, Erlykin and Wolfendale [7j. In this context it has recently been 
proposed by Salvati and Sacco [IS] that what MILAGRO sees is a signature of 
the supernova explosion that resulted in the formation of the Geminga pulsar. 
But one would expect all these processes to give rise to low-order multipole 
components, and indeed predominantly a dipole component, in the angular 
distribution. The peculiarity of the MILAGRO result is that it shows dis- 
crete features on small angular scales, and indeed not just random power in 
high-order multipoles, but quite discrete localised features. The data analysis 
technique used by MILAGRO does correspond to a high-pass filter which em- 
phasises these local hot-spots, but they cannot be dismissed as mere artifacts 
of the analysis. 



2 Neutron source models 



Motivated largely by the coincidence between the heliotail direction and region 
A, we first consider the possibility that region A might be explicable in terms 
of secondary neutron production in the gravitationally focused tail of ISM 
material that forms downstream of the Sun as it moves through the local 
ISM. Unfortunately, while the target has exactly the right size and location, 
the density is much too low to account for the observed signal (quite apart 
from the fact that this model would fail to account for region B even were it 
to work for region A). 



2.1 Neutron induced enhancement 



Observing into a narrow cone of solid angle dQ, a thin gas target of number 
density at distance r bathed in an isotropic and uniform flux of cosmic 
rays of number density ucr will produce a flux of neutrons 

r 2 dr dVt n T a pN cn CR 
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where <t p n is the cross-section for converting a proton to a neutron, either 
via the delta resonance or by directly knocking a neutron out of an atomic 
nucleus. The same flux of cosmic rays will produce a background rate within 
the cone of 

-^cn CR . (2) 
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Thus the relative enhancement of the signal will be simply 
/ n T a pN dr 



(3) 



and the question is whether we can get enough target material into the cone 
to give an enhancement factor of 10~ 3 . The ratio of the flux of secondary 
neutrons to cosmic rays can be estimated as [5] 



where the first factor is the so-called spectrum weighted moment for neutron 
production, o~ pp is the inelastic proton-proton cross-section and 77 ~ 1.5 is a 
correction factor to account for the presence of heavy nuclei. We note that 
while the production energy spectrum of the neutrons will be almost exactly 
the same as that of the ambient cosmic rays, the energy dependent (oc E) 
decay time will distort the spectrum and make it harder at certain energies, if 
the source is located at the "right" distance from us. For the estimated energy 
of the shower initiated particles of ~ 10 TeV, the source of neutrons should not 
be much further (otherwise the flux would be dramatically suppressed) and not 
much closer (otherwise the effect of spectral hardening would disappear) than 
d = tdecay(10 TeV)c = 7 n t n c ~ 0.1 pc where 7 n is the Lorentz time-dilation 
factor of the neutron and t n is the free neutron decay time. 



2.2 Gravitational focusing of ISM material 

It has long been known that the gravitational field of the sun can focus inter- 
stellar neutral atoms (and dust grains) as they flow past the Sun leading to 
the formation of a density caustic on the downstream symmetry axis. The ex- 
istence of such an effect has been directly confirmed by observations with the 
Ulysses spacecraft and indirectly by optical back-scatter observations [20~fT7] . 
Detailed analytical and numerical studies of this effect were first published by 
Fahr and co-workers back in the early 1970s [EH]. 

For our purposes a simple impulse approximation is adequate. A neutral par- 
ticle moving past the Sun at impact parameter b and velocity v experiences an 
acceleration of order GM/b 2 lasting for a time of order 2b/ v. It thus acquired 
a radial inward drift velocity of order 




(4) 



u = 



2b GM 
v b 2 



2GM 
vb 
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and reaches the symmetry axis after a time b/u — vb 2 /2GM and at a down- 
stream distance of 



L 



2GM 



(6) 



Clearly in the approximation of cold particles the density becomes infinite on 
the axis with a sharp cusp, but in reality the effect will be blurred by the 
initial velocity dispersion of the particles and in any case we can only observe 
with finite angular resolution. 

Let us therefore consider observing into a narrow cone of opening angle $ and 
consider first the case where the gravitational deflection is greater than At 
a focusing distance downstream of L the cone has a radius of dL and the flux 
of deflected particles entering the approximately cylindrical section between L 
and L + dL is nov2nbdb (where b and L are related as above) and the particles 
spend a time of order 2$L/u in the cone. The mean number density in the 
cone at distance L is thus 



for d <b/L. Once the deflection becomes comparable to the opening angle of 
the cone this approximation breaks down; there is still a focused component on 
axis, but it becomes less and less important relative to the directly inflowing 
material and the mean density in the cone tends asymptotically back to the 
ambient value n . That the enhancement factor is about 4 when the deflection 
angle equals the cone opening angle is intuitively obvious because in this case 
the collecting radius is precisely twice the radius of the circular cross-section 
of the cone into which the particles are forced. 

The peak column density is thus of order 4uqL where L is determined by 
the condition that b ~ dL. Taking d « 10~ 2 for an opening angle of about a 
degree and substituting values appropriate for the sun [13J (i;~3x 10 6 cms -1 , 
2GM/c 2 = 3 km) we get a maximum value of b of 3 x 10 15 cm and thus 
L = 3 x 10 17 cm = 0.1 pc. In fact the focusing will be probably be disrupted 
by the gravitational field of other near-by stars (at distances of a few parsecs) 
on scales any larger than this. 

2.3 Flux estimates 

We thus rather naturally produce a thin gas target at the appropriate distance 
of 0.1 pc, which can fit the location of region A and the harder spectrum, but 
unfortunately the effective enhancement of the ambient density by only a 
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factor of four is far too weak to account for the observed signal, the density of 
the local interstellar medium around the sun being quite low [13] with a rather 
accurately measured Helium number density of 0.015 cm -3 corresponding to 
uq ~ 0.15 cm -3 and thus a column density of 1.8 x 10 17 cm -2 giving a mere 
4 x 10~ 10 enhancement rather than the 10~ 3 observed. 

More generally we can conclude that no local neutron source model based 
purely on secondary production is viable because the required target density 
is so high, at some 10 6 times the local density that it should have been easily 
seen by astronomers studying the local bubble. 



3 Proton and ion source models 

If not neutrons, let us now consider the possibility that the signal might be 
due to protons (or more generally heavy ions). Even supposing the Earth 
to be magnetically connected to a remote accelerator, and for ions to travel 
rather freely along this field, one would naively expect them to arrive with 
a broad distribution of pitch angles relative to the field and it is thus hard 
to explain the narrow angular scale of region A (the angular distribution of 
arrival directions is quite compact with a size of only a few degrees, whereas a 
gyrating streaming proton distribution should fill virtually the entire forward 
hemisphere of 2ir steradians). 

This objection applies a fortiori to the model in [18] where the interstellar 
propagation is assumed to be strong-scattering diffusion in the Bohm limit. On 
this model the anisotropy signal would be of dipole form with significant back- 
scattering of particles. In fact, while Bohm scaling for the scattering in the 
immediate vicinity of the shock is quite plausible and indeed seems required, 
it would be extraordinary for interstellar propagation to be described in this 
way and it would completely contradict the conventional picture of interstellar 
magnetic fields and cosmic ray propagation. 

Similar considerations apply to any local accelerator of protons using stochas- 
tic acceleration (first or second order fermi) all of which produce essentially 
isotropic particle distributions. Magnetic reconnection can produce tightly col- 
limated beams of particles along the magnetic null line, but it is unlikely that 
the fields can be so well organised over the scales required in the solar neigh- 
bourhood. In any case there is the fundamental problem that any such local 
acceleration process is constrained to the Hillas limit on attainable particle 
rigidity of BVL where B is the magnetic field strength, V a characteristic ve- 
locity and L a characteristic length. Taking optimistic values of a 2 /iG field, 
a length scale of 0.1 pc and a velocity of 30kms _1 (the order of magnitude of 
the solar velocity with respect to the LISM) gives a maximum particle rigidity 



5 



of 20 GV, well short of the energies at which the signal is seen. And of course if 
the source is less than a gyroradius away, by definition the accelerator cannot 
use magnetic confinement anyway and we would have to find an electrostatic 
acceleration potential of 10 13 V in the local neighbourhood! In this context we 
note that Salvati and Sacco [H] also give an energy budget argument against 
a local heliospheric origin which we find compelling. 

Because the essential problem is the small gyroradius one could try using 
particles with a larger mass to charge ratio. In theory one could consider 
singly charged heavy ions (similar to the anomalous cosmic rays observed 
at much lower energies), but even if heavy ions could be accelerated while 
retaining an electron cloud, the electron stripping cross-sections are so large 
at these energies that any such particles would almost instantaneously become 
fully ionised during propagation. In any case, even assuming the signal to be 
singly charged heavy ions would only gain a factor of one to two decades and 
thus still leave the putative accelerator closer than the nearest stars. Charged 
dust grains can only be slightly accelerated before frictional and sputtering 
lossses intervene jH] and in any case on interaction with the atmosphere are 
not expected to produce the penetrating showers seen by MILAGRO. 

One possible way out would be to suppose that we happened to be connected 
to a transient event and only the fastest particles (moving essentially straight 
along the field with little or no gyration) have reached us, but this seems very 
contrived (although it does make the interesting prediction that the spot will 
brighten and broaden out on a time scale of the light travel time to the event 
which could be as short as a few years). 

A more plausible alternative is to invoke relatively scatter-free propagation 
with magnetic moment conservation, and to suppose that the field line con- 
necting us to the source has come from a region of high field into a much 
lower field region (a rather plausible scenario). Adiabatic conservation of the 
magnetic moment then implies that the particles momentum component per- 
pendicular to the field also has to drop in such a way that p±/B is constant 
and thus the particles can be focused into a beam parallel to the field. An- 
other way to think of this is in terms of Liouville's theorem on the conservation 
of phase-space density. To compensate for the reduction in angular distribu- 
tion from 2-7T steradians down to the observed 0.3 steradians of regions A and 
B, the cross-section of the magnetic flux tube has to expand by an amount 
27r/0.3 ~ 20 which seems quite plausible. 

Of course if there is a significant flux of such particles this is a highly unstable 
velocity distribution and resonant modes will be excited which scatter the 
particles and attempt to restore isotropy, but for a sufficiently weak beam 
this effect will not be a problem and the resonant modes can be damped 
by, among other possibilities, Landau damping on the bulk cosmic rays and 
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ion-neutral friction [T91ITT] (and we are only talking about a beam of order 
1CT 3 relative to the background). In addition, if the field lines continue to 
diverge, the adiabatic focussing effect can at least partially compensate for 
pitch-angle scattering as long as the characteristic length of the field variation 
is comparable to the scattering mean free path. On balance this seems the 
least unlikely explanation of the results. A field configuration departing from 
strict axisymmetry could allow the loss-cone to be defocused into the two 
regions observed by MILAGRO, but the alignment with the heliotail of region 
A would have to be largely coincidental on this explanation. However it is 
possible that draping of the interstellar field along the heliosheath might help 
and this may also account for the fact that, as noted by Abdo et al pQ, the 
directions of regions A and B do not seem to lie along what is thought to be 
the local interstellar field direction. But it is not clear that the local field is 
really know at the scales of interest here; to quote Leroy [T2] "The interstellar 
magnetic field revealed by the nearest polarized stars is certainly not a smooth, 
uniform field along the Galactic equator" . 



4 Conclusion 

The MILAGRO observations of localised small-scale structures in the arrival 
direction distributions of cosmic rays represent a real challenge to our under- 
standing of cosmic ray origin and propagation. It would appear that the only 
viable explanation involves relatively scatter-free focused transport along a 
diverging magnetic field connecting the solar neighbourhood to a magnetic 
mirror where the field increases by about a factor of 20 and where the cosmic 
ray intensity on the other side of the mirror is slightly higher in amplitude 
and harder in spectrum. The particles in regions A (and presumably B?) are 
then those particles in the so-called "loss cone" which penetrate through the 
trap. The bulk of the cosmic rays we see have been mirrored off the trap 
and reflect the general population in the solar neighbourhood. This is illus- 
trated in cartoon form in Fig. [TJ We note in passing that relative motion 
between the sun and the trap (presumably a field concentration associated 
with a nearby molecular cloud) could also induce a large-scale signal in the 
anisotropy data. Further work, which must also address the question of the 
large scale anisotropy seen by MILAGRO as well as by the Tibet Air Shower 
Array and the Super-Kamiokand experiments is clearly required to un- 
derstand these fascinating results. 

Finally, we note that while this model allows the small angular scale features 
to be due to a source located much further away than one proton gyroradius, 
it is still advantageous to have the source relatively near-by, say at less than 
100 pc. In fact as estimated by a number of authors, such a single nearby source 
with a total energy release of about 10 50 erg can significantly contribute to the 
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Fig. 1. A cartoon representation of the suggested model. The magnetic field near 
the sun is linked through a strong field region where most particles are mirrored 
to a region containing a slightly higher amplitude and harder spectrum of cosmic 
rays, presumably as a result of a local source event such as a supernova explosion. 
Only particles with pitch angles closely aligned to the field can penetrate through 
the mirror and emerge as a field-aligned beam on the other side (another possibility 
would be to locate the source inside the strong field region itself rather than behind 
the trap). 

observed cosmic ray flux, see for example references [2ir7fT0] , and in fact there is 
no doubt that there must have been at least one recent nearby cosmic TeVatron 
from the detection of cosmic ray electrons at E > ITeV [15J. Interestingly a 
detectable cosmic ray anisotropy from such a source has been predicted [TU] 
and it is likely that the LAT on GLAST will be able to provide detailed studies 
of such TeV electrons [U]. Without wishing to claim any particular event as 
the source, the suggestion that the supernova explosion which resulted in the 
formation of the Geminga pulsar may have been sufficiently close and recent 
to leave detectable traces we find quite plausible and in this we agree with 
Salvati and Sacco [IB]. Other candidates, as pointed out by Kobayashi et al 
[TO] who discuss this extensively, are the Vela, Cygnus loop and Monogem 
remnants. 
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